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To study the nature of partially substituted Yb-ions in a Ce-based Kondo lattice, we fabricated 
high quality Cei-^Ybj; Coins epitaxial thin films using molecular beam epitaxy. We find that the 
Yb-substitution leads to a linear decrease of the unit cell volume, indicating that Yb-ions are divalent 
forming Kondo-holes in Cei-aYb^CoIns, and leads to a strong suppression of the superconductivity 
and Kondo coherence. These results, combined with the measurements of Hall effect, indicate 
that Yb-ions act as nonmagnetic impurity scatters in the coherent Kondo lattice without serious 
suppression of the antiferromagnetic fluctuations. These are in stark contrast to previous studies 
performed using bulk single crystals, which claim the importance of valence fluctuations of Yb-ions. 
The present work also highlights the suitability of epitaxial films in the study of the impurity effect 
on the Kondo lattice. 

PACS numbers: 72.15.Qm, 74.62. Dh, 74.70. Tx, 81. 15. Hi 



The impurity effect has been one of the key issues for 
the study of unconventional superconductors including 
cuprates EL iron-pnictides Q and heavy fermion com- 
pounds [3|— lOf] . In these strongly correlated electron sys- 
tems, the superconducting state emerges from the com- 
petition between different phases, such as a magnetically 
ordered phase. Understanding the impurity effect is im- 
portant because it is intimately related to the supercon- 
ducting pairing interaction and competing electronic cor- 
relations. In heavy fermion compounds, it was generally 
found that a substitution of magnetic sites with nonmag- 
netic impurities yields a substantial reduction in both the 
onset of coherence temperature T co h at which the for- 
mation of heavy fermions occurs, and the superconduct- 
ing transition temperature T c . Recently, however, highly 
unusual impurity effect has been reported in the heavy 
fermion compound CeCoIns, which is in sharp contrast 
to the above tendency [6H8J]. 

CeCoIns with tetragonal symmetry EH is a very clean 
superconductor with large mean free path, which is 
known to exhibit a plethora of fascinating electronic 
properties. The unconventional superconductivity with 
d-wave symmetry app ears in the vicinity of magnetic 
ordered phase Ill4l3j. The normal state exhibits pro- 
nounced non-Fcrmi-liquid behaviors, which are believed 
to be due to the proximity of an antiferromagnetic (AF) 
quantum critical point [14j, |_15 1 . Recent study revealed 



that the superconductivity occurs in two-dimensional 
Kondo lattice composed of square lattice of Ce atoms 
EH- Thus CeCoIn 5 may provide an ideal playground 
for investigating the impurity effect on unconventional 
superconductivity in strongly correlated /-electron sys- 
tems. The impurity effect on CeCoIns has been inves- 
tigated extensively by rare earth R (or Hg, Cd, Sn) 
substitution for Ce (or In) [!, 0, El El- Two notable 



features have been pointed out. Firstly, non-magnetic 
impurities locally suppress superconductivity, generat- 
ing an inhomogeneous electronic "Swiss cheese" , simi- 
lar to cuprates [B|. Secondly, several groups have re- 
ported more striking and unexpected results in the case 
of Yb-substitution on Ce-site: the superconductivity in 
Cei_ a; Yb a ;CoIn5 is robust against Yb-substitution and 
T c decreases linearly with x towards K as x — > 1, while 
T co h remains essentially unaffected with Yb-doping 0- 
8]. YbCoIns is a conventional nonmagnetic metal with 
no superconducting transition down to 20 mK, indicat- 
ing that Yb is close to divalent and in the nonmagnetic 
closed shell 4/ 14 configuration 1(1, T^. Therefore, the ro- 
bustness of the superconductivity upon Yb-substitution 
in CeCoIns with unconventional pairing symmetryis ex- 
traordinary. Moreover, according to Refs.Q and El, the 
lattice parameters remain roughly constant as x changes 
until the phase separation occurs at x ~ 0.8, indicating 
the violation of Vegard's law. These results are markedly 
different from other R substituted CeCoIns, in which R 
substitutions suppress the superconductivity at approxi- 
mately x = 0.2-0.3 E3, Ull- These anomalous behaviors 
have been discussed in the light of the valence fluctua- 
tions of Yb-ions However, the electronic state of Yb 
in the Kondo lattice of CeCoIns is poorly understood. 

To obtain further insight into the nature of Yb- 
ion in the Ce-based Kondo lattice, we fabricated 
Cei-^YbzCoIns thin films using molecular beam epi- 
taxy (MBE) technique. This technique, recently ad- 
vanced by our group, enables the growth of uniform 
heavy fermion thin films, as demonstrated by the suc- 
cessful growths of Celn.s /Laln^ 21| and superconduct- 



ing CeCoIn 5 /YbCoIn 5 [3| superlattices with one- unit- 
cell thick Celn3 and CeCoIns layers, respectively. Since 
the epitaxial growth occurs in non equilibrium condition 
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FIG. 1. (color online), (a) X-ray diffraction 0-scan data of 
the (115) peak for Ceo.sYbo^CoIns. Inset : Streak patterns 
of the RHEED image during the crystal growth, (b) X-ray 
reciprocal lattice mapping for Ceo.sYbo^CoIns near the (115) 
peak, (c) Lattice parameters a (filled green triangles) and c 
(filled red squares and circles) for Cei_a; Yb^CoIns films. Here 
c is determined by reciprocal mapping (squares) and 9 — 29 
scan (circles). Open blue circles represent c of bulk crystals 
reported in Ref. [Q]. The dased lines are guides for eyes. The 
inset shows x determined by EDX, iedx, vs. x determined 
by the evaporation ratio of Ce to Yb atoms. 



at temperatures much lower than the temperature used 
for single crystal growth by flux method, MBE is suitable 
for the preparation of systems with homogeneously dis- 
tributed Ce and Yb-atom. Moreover, owing to the ability 
to evaporate Ce and Yb atoms simultaneously, the ratio 
of Ce to Yb can be controlled precisely. Using our thin 
films, we find that the Yb-substitution leads to a strong 
suppression of T c , in accordance with Abrikosov-Gorkov 
(AG) theory. These results are in sharp contrast to pre- 
vious studies performed using bulk single crystals [H-df. 

The c axis oriented epitaxial Cei-^Yb^CoIns films 
were grown by MBE. The (001) surface of MgF 2 with 
rutile structure (a = 0.462 nm, c = 0.305 nm) was used 
as a substrate. The substrate temperature was kept 
at 530-550 °C, depending on the composition. Each 
metal element was evaporated from individually con- 
trolled Knudscn-ccll. The typical deposition rate was 
0.01-0.02 nm/s. The area and thickness of the films were 
5.0 x 10.0 mm 2 and 120 nm, respectively. The epitax- 
ial growth of each layer with atomic flatness was care- 
fully checked by monitoring the streak patterns of the 
Reflection High Energy Electron Diffraction (RHEED) 
during the deposition (the inset of Fig. la). The X- 
ray diffraction (/)-scan with 4-fold peaks shown in Fig. la 
also indicates the epitaxial growth of the film, which was 
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FIG. 2. (color online). Temperature dependence of the (a) 
resistivity and (b) Hall coefficient for Cei-^YbzCoIns. 



further confirmed by the X-ray reciprocal lattice map- 
ping as shown in Fig lb. Surface flatness in atomic level 
was confirmed by the atomic force microscopy. Excellent 
chemical homogeneity across the whole sample area was 
confirmed by the energy dispersive X-ray spectroscopy 
(EDX) analysis. We could not grow the epitaxial films 
for x > 0.4 without buffer layers. To grow YbCoIns epi- 
taxial films, we used Celn3 as buffer layers on the sub- 
strate. As shown in the inset of Fig. lc, the ratio of Ce to 
Yb of the measured films determined with EDX coincides 
well with the evaporation ratio of Ce to Yb atoms. 

In our epitaxial thin films (a = 0.462 nm, 
c = 0.753 nm), lattice parameter a is dictated by the sub- 
strate lattice parameter and misfit strain slightly enlarges 
a from that of bulk single crystal value (a = 0.461 nm, 
c = 0.755 nm). Figure lc shows the doping evolution of 
lattice parameters determined by X-ray diffraction. Be- 
cause of the epitaxial strain effect, a is independent of 
x and coincides well with a of MgF 2 substrate. On the 
other hand, c (and hence unit cell volume) decreases lin- 
early with increasing x and lies on the line connecting 
x = and 1. According to the Vegard's law, the unit cell 
volume should decreases linearly with Yb-substitution, if 
there are no changes in the valence of Ce and Yb ions. 
Therefore the observed linear relation implies that the Ce 
and Yb ions in the present Cei-^Yb^ Coins films retain 
the valences of the end member compounds, i.e. the Yb- 
ion is divalent in Cei-^Yb^CoIn.j, forming nonmagnetic 
"Kondo holes" in the 4/ 1 lattice. This result is in sharp 
contrast to the results using bulk single crystals, in which 
c is x- independent (open circles in Fig. lc) 

Figure 2a and 2b depict the temperature dependence 
of the in-plane resistivity p and Hall coefficient i?u for 
H\\ c, which is defined as the field derivative of Hall 
resistivity at the zero field limit, respectively. First we 
compare the transport properties of pure CcCoIns thin 
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FIG. 3. (color online), (a) x-dependence of T co h normalized 
by the value at x=0 for thin films (filled red circles) and bulk 
crystals reported in Ref. 0] (open blue circles), (b) Tempera- 
ture dependence of the resistivity at low temperatures. The 
dased lines are the extrapolation above T c . 



film with those of single crystal. The temperature depen- 
dence of both the resistivity and Hall coefficient in thin 
film is essentially the same as those in bulk single crystals 
[lit . Associated with an incoherent-coherent crossover, 
p shows a maximum at Tcoh at around 40 K. The abso- 
lute value of p of CeColns thin film is close to that of 
bulk single crystal. In particular, p ~ 6 p,Q cm at the 
superconducting onset is close to the bulk single crys- 
tal value, which ranges from 3 to 8 pfl cm [3I I22L [23j|. 
The resistivity below T co h in bulk single crystal exhibits 
a T-linear dependence, which is a hallmark of non-Fermi 
liquid behavior. In our thin film, the resistivity exhibits 
a T-dcpendence with exponent slightly below unity. 

i?H is negative in the whole temperature range. The 
temperature dependence of Rh is closely correlated with 
that of the resistivity [TH, I24J . At high temperatures 
T > T co h, Rh is nearly temperature independent. Below 
Tcoh, -Rh decreases rapidly as the temperature is low- 
ered. Further decrease of temperature increases i?n after 
showing a minimum at around 5 K. The magnitude of i?H 
above T co h matches well with that of bulk single crystal 
[lij . The enhancement of the absolute value of i?H be- 
low T co h is nearly 15 times that at above T co h, which 
is nearly half of that in bulk single crystal. The super- 
conducting transition temperature of the CeCoIn 5 thin 
film is 1.95 K, which is slightly lower than T c = 2.3 K 
of bulk single crystal. Since the p value at the super- 
conducting onset (and residual resistivity which will be 
discussed later) in thin film is comparable to the bulk 
single crystal value, this T c reduction is likely to be due 
to the strain effect arising from the substrate. In fact, it 
has been reported that the negative pressure in CeCoIn 5 
reduces T c 25[. Based on these results, we conclude that 
the quality of our epitaxial thin films grown by MBE is 
comparable to that of high quality bulk single crystals. 

Next we discuss the effect of Yb-substitution on the 



normal state transport properties. As shown in Fig. 2a, a 
broad maximum of p(T) associated with the formation of 
Kondo coherence is observed at T co h(x) in Yb-substituted 
CcCoIn5. Figure 3a depicts the x-dependence of T co h, 
which is simply defined as the peak position of p(T). Yb- 
substitution seriously reduces the Kondo coherence and 
T C oh appears to go to zero at around x = 0.5. Since it 
is highly unlikely that Fermi temperature is dramatically 
suppressed by Yb-substitution, the observed reduction of 
T C oh arises from the destruction of Kondo coherence by 
the Yb 2+ Kondo- holes in the Ce 3+ /-electron lattice. In 
Yb-dopcd compounds p(T) exhibit sub-T-linear behavior 
below T co h- We try to fit the resistivity by a power law 
p(T) = pa + AT a (Fig. 3b). The resistivity can be fitted 
very well by this equation for x — and 0.02 with the 
residual resistivity po= 4.8 and 6.0 /zf2 cm and a=0.88 
and 0.80, respectively, in a wide temperature range below 
~ T coh /2 down to T c . For x > 0.05, we cannot fit p(T) 
by this equation in a wide temperature range and we 
then estimate po by the polynomial fitting above T c to 
Tcoh/2 (see the dashed lines in Fig 3b). As shown in 
Fig. 4a, the residual resistivity nearly linearly increases 
with x. The present results are very different from those 
reported in bulk single crystals, in which T co h is nearly 
x-independent as shown in open symbols in Fig. 3a and 
a barely changes with x up to x — 0.1 jf|. 

The measurements of Hall effect provide vital informa- 
tion for the nature of the Yb-ion in the Ce-Kondo lattice. 
As shown in Fig. 2b, the Hall effect is very sensitive to 
the Yb-substitution. At T > T co h, the absolute value of 
i?H is small and nearly temperature independent for all 
x. This is consistent with the the Vegard's law which 
suggests divalent Yb in trivalent Ce lattice. Below T co h, 
the enhancement of |-/?h| is strongly suppressed by the 
Yb-doping. The low-temperature enhancement of |-/?h|, 
which is significantly larger than [1/riej where n is the 
carrier concentration, have also been reported in other 
strongly correlated electron systems including cuprates 
[H iron-pnictides |27|, K-(BEDT-TTF) 2 Cu(NCS) 2 [HI 



and V2- y 03 [29j. It has been shown that the striking en- 
hancement of |i?n| at low temperatures can be accounted 
for in terms of "hot spot" on the Fermi surface (FS) and 
backflow effect, both of which originate from strong AF 
fluctuations [30J]. The hot spot is a FS region where the 
electron lifetime is unusually short, which appears at the 
positions where AF Brillouin zone boundary intersects 
with the FS. Since the hot spot docs not contribute to 
the electron transport, the effective carrier concentration 
is reduced, which results in the enhancement of [-Rh| from 
|l/ne| 3l|. However, this effect is not enough to explain 
the dramatic enhancement of |-Rh| below T co h- Another 
important effect is the backflow effect accompanied by 
the anisotropic scattering, which is called the current 
vertex correction [3(3]. Backflow is a polarized current 
caused by a quasiparticle excitation inherent in Fermi 
liquid. In the presence of the backflow effect, the total 
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FIG. 4. (color online), (a) s-dependence of the residual re- 
sistivity. Dotted line is a guide for eyes, (b) rr-dependence 
of the superconducting transition temperature normalized by 
the value at x=0 for thin films (filled red circles) and bulk 
crystals (open blue squares @] and open circles Q|). (c) T c /T c o 
plotted as a function of residual resistivity. Dotted line is the 
linear extrapolation from the low po region. The solid curve 
represents the result fitted by Abrikosov-Gorkov (AG) theory. 



current is not parallel to the Fermi velocity and then the 
Hall coefficient is seriously modified from the Boltzmann 
value derived from the curvature of the FS. It has been 
shown that the backflow effect can largely enhance \Ra\ 
in CeMIn 5 (M=Co, Rh and Ir) [H[24|. 

The observed strong suppression of \Rn\ by Yb- 
substitution indicates that the backflow effect is seriously 
reduced. There are two possible origins for this. The 
first is the reduction of AF fluctuations and the second 
is the increase of isotropic impurity scattering by non- 
magnetic Yb-ions. In the former case, since the system 
is tuned away from the quantum critical point, the T- 
linear resistivity is expected to change to T 2 -dependence 



accompanied by a reduction of po [15| . However, such a 
trend is not observed here. Instead, the T-linear resis- 
tivity changes to sub- T-linear and po increases with Yb- 
substitution. This increase of po is consistent with the 
increase of isotropic scattering. In addition, it has been 
shown that in the presence of strong AF fluctuations, 
sub-T-linear dependence of p(T) appears with increas- 
ing impurity scattering [32*] . which is again consistent 
with the present results. Thus the normal state trans- 
port properties suggest that Yb-ions are divalent and act 
as impurity scattering centers, rather than a suppressor 
of AF fluctuations. 

Next we discuss the influence of Yb-substitution on the 
superconducting properties. As shown in Fig. 3b, sharp 
resistive transitions are observed for x < 0.20, but the 
transition is slightly broader for x = 0.30. No super- 



conducting transition is observed for x = 0.40 down to 
0.3 K. Figure 4b displays the x-dependence of T c , which 
is defined as the midpoint of the resistive transition, nor- 
malized by T c of x=0. T c decreases linearly with x and 
goes to zero at around x = 0.40. We stress that the ob- 
served Yb-doping evolution of T c is again in marked con- 
trast to that reported in bulk single crystals shown by 
open symbols in Fig. 4b. In Fig. 4c we plot T c /T c q (see 
below) against po- In the low po range, T c /T c o decreases 
linearly as shown by the dotted line. When po exceeds 
~ 18 pfl cm, T c begins to deviate from linearity. We 
analyze this trend of T c /T c q in accordance with the AG 
theory of nonmagnetic impurity effect in ri-wave super- 
conductors [33| . The fitting parameters used to compare 
experiment to theory are the slope in the low p region 
and T c q, where T c q is the transition temperature with no 
pair-breaking. The solid line in Fig. 4c is the result of fit- 
ting obtained by using T c0 = 2.20 K and the initial slope 
-0.025 (ptt cm) -1 . As shown in Fig. 4c, the suppression 
of T c is well reproduced by the AG pair breaking curve 
in the whole po range. 

The fact that our data can be well described by the 
AG theory immediately implies that the Yb-ions are ran- 
domly distributed with weak inter-ion correlation. More- 
over, the "Swiss cheese" model [H| is inadequate because 
such a model will give a T c (p )-curve that lies above the 
dotted line for high p in Fig. 4c |34( . The effect of Yb- 
substitution on CeCoIns bears striking resemblance to 
other rare earth substitutions, although Yb is divalent 
while other rare earths are trivalent. This suggests that 
these ions act as impurity centers with unitary scatter- 
ing, regardless of their valence. Furthermore, the va- 
lence fluctuations of Yb-ions, if present at all, do not 
appear to play an important role on the physical prop- 
erties in the normal state as well as the superconducting 
state. At the present stage, the essential difference of the 
Yb-substitution effect between thin films and bulk sin- 
gle crystals remains an open question. A possible origin 
for this may be that bulk crystals contain some regions 
where inter- Yb- ion correlation is important (35| . Further 
studies are required to clarify this issue. 

To summarize, in high quality Cci-^Ybj-CoIns epitax- 
ial thin films, Yb-ions are divalent and no signature of 
the valence fluctuation is observed. The Yb-substitution 
leads to a strong suppression of the superconductivity 
and Kondo coherence. The suppression of T c can be 
well described by AG theory. These results indicate that 
Kondo-holes created by Yb-ions act as nonmagnetic im- 
purity scatters in the Kondo lattice with no serious re- 
duction of AF fluctuations. These are in sharp contrast 
to previous studies performed using bulk single crystals, 
which claim the importance of valence fluctuations of Yb- 
ions. The present work also emphasizes the uniqueness 
of the epitaxial films in the study of the impurity effect 
on the /-electron Kondo lattices, due to the prospect of 
preparing highly homogeneous doped systems. 
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